Granular spirals on erodible sand bed submitted to a circular fluid motion 
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An experimental study of a granular surface submitted to a circular fluid motion is presented. The 
appearance of an instability along the sand-water interface is observed beyond a critical radius Vc- 
This creates ripples with a spiral shape on the granular surface. A phase diagram of such patterns 
is constructed and discussed as a function of the rotation speed lu of the flow and as a function of 
the height of water h above the surface. The study of Vc as a function of h, uj and r parameters is 
reported. Thereafter, rc is shown to depend on the rotation speed according to a power law. The 
ripple wavelength is found to decrease when the rotation speed increases and is proportional to the 
radial distance r. The azimuthal angle e of the spiral arms is studied. It is found that e scales 
with huir. This lead to the conclusion that e depends on the fluid momentum. Comparison with 
experiments performed with fluids allows us to state that the spiral patterns are not the signature 
of an instability of the boundary layer. 

PACS numbers; 45.70.Mg, 81.05.Rm, 64.75.+g 



I. INTRODUCTION 

Many experiments have been dedicated to the study 
of fluid flows between rotating stainless-steel disks |jl|, 
^, Q . Three experimental conflgurations can be found in 
the hterature: two contra-rotative disks; the upper disk 
rotating while the other is at rest; and the lower disk ro- 
tating, the other being at rest. For a given Reynolds 
number, i.e. at a given distance from the center of 
the plate a transition in the boundary layer occurs. 
The encountered instabilities are of Karman, Ekman or 
Bodewadt types and may lead to one of the following pat- 
terns: (i) concentric circles at a low rotation speed (ii) 
spirals at a medium rotation speed and (iii) a disordered 
state at large rotation speeds Q . 

On the other hand, one can note the growing interest of 
the statistical physics community to the granular state. 
Particularly, the ripple formation on sand bed eroded by 
a fluid such as air or water has become a largely studied 
phenomenon |7[ ||, |lO| . Despite the familiar aspect 
of the ripples, the physical mechanisms involved in their 
formation are related to complex phenomena of granu- 
lar transport such as avalanches, saltation, reptation and 
suspension [Tl| . 

In the present paper, both subjects will be combined 
in a single experiment. We will study the stability of a 
granular surface submitted to a circular fluid flow. In the 
next section, the experimental setup will be described. In 
Section III, we will present the results, and discuss them 
in Section IV. Finally, a summary of our findings is given 
in Section V. 



II. EXPERIMENTAL SETUP 

Our experimental setup is illustrated in Fig.|l| and con- 
sists of a horizontal circular plate connected by a belt 
to an engine. The rotation speed can be adjusted from 
8 rpm up to 100 rpm. A cylindrical container (radius 



i? = 11 cm) filled with water and sand (mean grain diam- 
eter d = 282.5 //m) is placed in the center of the plate and 
put into rotation. When the rotation is brutally stopped, 
the sand bed remains fixed while the water continues its 
inertial circular motion. The shear stress applied to the 
sand by the water flow initiates various grain motions: 
saltation, reptation and suspension. After a short time 
(typically 2s), current ripples are formed [see Fig.0 bot- 
tom]. 
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FIG. 1: (Top) Experimental setup for the ripple formation 
and analysis. (Bottom) Profile of the landscape along the cir- 
cumference of the container after the formation of the ripples. 
Notice that the image has been extended vertically for clarity. 



A CCD camera is placed on the top of the plate and 
records an image of the sand bed. In order to avoid 
shadow effects, the container circumference is homoge- 
neously illuminated. An annular luminescent tube with 
a radius of 12 cm is placed horizontally around the con- 
tainer, at the level of the sand-water interface. We as- 
sume that the height of sand is given by the intensity of 
the pixels on the grayscale images, due to this side light- 
ning. Within this assumption, the landscape is given 
with an accuracy of 0.7 d « 200 fim. 

The parameters of our experiment are: (i) the height 
of water h over the sand bed before rotation, and (ii) the 
rotation speed lu of the container. Moreover, one should 
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note that the Reynolds number varies on the granular 
surface as a function of the distance < r < R from 
the center of the container H]. Therefore, the Reynolds 
number Re we will use can be written as 

R. = (1) 

where v is the kinematic viscosity of the water, i.e. v = 
10"^ m^ s-\ 



III. RESULTS 

We have performed experiments with values of lu rang- 
ing from 29 rpm to 71 rpm, and values of h ranging 
from 1.05 cm to 3.36 cm. For each experiment, an im- 
age of the landscape has been taken. On each image, 
we have extracted transversal profiles of the landscape 
at 50 different r values, i.e. different R^ values. Figure |^ 
presents typical data for two values of the Reynolds num- 
ber: Re = 10.4 (r=3 cm) and Re = 34.6 (r=10 cm). For 
small Reynolds number values, the profile is very noisy 
and no modulation is observed. The noise comes from 
the irregularities of the granular surface. On the con- 
trary, we have noticed a periodic profile for large values 
of Re- The transition between noisy and periodic regions 
is well defined and occurs at a given critical radius re- 
The profiles are periodic above re and noisy elsewhere. 
This kind of behavior is characteristic of the emergence 
of an instability in the sand-water interface. This is sim- 
ilar to what is observed in the case of an instability in 
the boundary layer of a fluid flowing over a stainless-steel 
disk [|. 




FIG. 2: Height of sand as a function of the position at a given 
Reynolds number: (Top) Re = 10.4, only fluctuations due to 
sand grains are visible. (Bottom) Re = 34.6, the instability is 
characterized by a periodic profile. The experimental param- 
eters are: ui — 45.18 rpm and h — 1.6 cm. Note that images 
have been rescaled in order to be compared. 

For all tested values of both parameters h and w, we 
have noticed the appearance of that instability for high 
Reynolds number values. However, we have remarked 
that the pattern formed on the granular landscape de- 
pends on the experimental parameters. Indeed, we have 
recorded three kinds of patterns : (i) one spiral; (ii) two 
superposing spirals, one having the arms close to the tan- 
gential direction to the container circumference, the other 



being more open; and (ii) two paired spirals having the 
same number of arms and having their arms coupled one 
to the other. Figures ^-c illustrate those different pat- 
terns. 




FIG. 3: Depending on the rotation speed, three patterns 
can be observed: one spiral, two superposing spirals and two 
paired spirals, (a) At low rotation speeds {uj < 35 rpm) the 
ripples form one spiral, (b) At medium speeds (35 < w < 45 
rpm) two spirals are superposed, one having arm azimuthal 
angles larger than the other ones, (c) At high rotation speeds 
{u > 45 rpm) the arms of the two spirals are coupled. 



We have constructed a phase diagram giving the pat- 
tern as a function of both to and h parameters. One 
should note [see Figure |^ how the pattern depends on 
the rotation speed lo and is less dependent on the height 
of water h. Looking for details in images corresponding 
to the phase diagram, we have concluded that the pat- 
terns evolve as follows. At a low rotation speed {lu < 35 
rpm), spirals are created above a critical radius Ve- An 
increase of lu causes the increase of the spiral arm an- 
gles, as we will see further. At medium rotation speeds 
(35 < w < 45 rpm) a second spiral appears. The arms 
of this second spiral are nearly tangential to the cylinder 
circumference. As a consequence, one observes two su- 
perposed spirals. If the rotation speed is increased again 
{lu > 45 rpm) , the second spiral becomes more and more 
open, so that both spirals are paired. 



A. Fourier analysis 

In order to determine the geometrical properties of the 
patterns, we have decomposed the data of the different 
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FIG. 4: Phase diagram giving the observed pattern as a function of the rotation speed uj and the height of water h. The rotation 
speed of the container clearly plays a role in the pattern formation, while the height of water has no significant effect. The 
three phases are separated by solid lines. At a; < 35 rpm, only one spiral is observed. A second spiral appears for 35 < tj < 45 
rpm and is paired to the former spiral for a; > 45 rpm. 



profiles within Fourier series, according to the equation 
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where f{x) is the height of sand at position x along the 
circle of radius r, m is the Fourier mode, a™ and bm 
are the Fourier coefficients. The sum is extended over 
A'^ terms, where TV is equals to the number of points in 
the discretization of f{x). Typically, N « 2000. One 
should note that the geometry of our experimental setup 
allows us to use Fourier approximation without restric- 
tion. From this approximation, we can compute the sta- 
tistical weight Am of a given mode m in the series by 



(3) 



The knowledge of all weights Am allows us to determine 
the mean amplitude (A) of the ripples 
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The mean wavelength of the patterns at the radius r is 
then given by 
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(A) = 27rr 



(5) 



which is the circumference at the considered radius r, 
divided by the statistical mean of the weighted modes 
m. It is important to note that practically, the sums 
are calculated from m — A instead of m = 1. In such 
a way, we remove the large modulations of the surface. 
Those large oscillations are artifacts due to the prepa- 
ration of the sand bed before the experiment is started. 
Such problem is also encountered in other experiments 
dealing with granular surfaces ^. 



B. Critical radius rc 

With the help of Fourier decompositions, one can eas- 
ily find the value of r^. For all radii r < r^, the profiles 
are noisy. In such a noise, all frequencies have the same 
probability of being found in the profiles (white noise). 
Mathematically, this fact results in Fourier series which 
are sums of sinus and cosinus with nearly the same sta- 
tistical weights Am- On the contrary, for r > Tc the 
periodic shape of the profiles implies the existence of a 
characteristic mode M in each profile. The value of M is 
equals to the number of ripples in the profile. The plot of 
Am as a function of the Fourier modes will thus exhibits 
a peak corresponding to the mode M. Eventually, the 
critical radius Vc is equal to the radius where the first 
peak is recorded. Practically, Tc is defined as the radius 
where the highest peak reaches a defined threshold value. 

In Figure ^, we report the amplitudes Am of the dif- 
ferent Fourier modes as a function of the radius r. One 
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can see that peaks appear only for r > 4.1 cm, meaning 
r„ = 4.1 cm. Notice that in this case, M ^ 6 ioi r — rr- 



22 




FIG. 5: Amplitude Am of the different Fourier modes m as a 
function of the radius r. Since the instability is characterized 
by periodic profiles, Fourier peaks are observed in the spec- 
trum only for radii r > Vc- On the contrary, stable regions 
have no periodic characteristics. The critical radius rc = 4.1 
cm is emphasized by a solid line. The experimental parame- 
ters are: uj = 45.18 rpm and h — 1 cm. 



We have studied the dependence of rc on the height 
of water and the rotation speed. Figure ^ presents the 
mean ripple amplitude (A) as a function of the radius r 
for different values of h at fixed oj — 41.38 rpm. First, 
one should note that the r-axis can be decomposed in 2 
distinct regions: (i) If r < Tc, all curves are nearly con- 
stant. This is the stable part of the sand-water interface, 
(ii) For r > Tc, all curves grow. Notice that the critical 
radius is the same for all curves. This means that the 
birth of the instability is not controlled by the height of 
water. Moreover, one can see how (A) is inversely pro- 
portional to h. Thus, the higher the velocity gradient of 
the water is, the larger are the ripples. 

Plotting the critical radius function of the ro- 

tation speed LU [see Figure |^ shows that the instability 
rises as sooner as the fluid speed is high. Fitting the data 
with the single power law 



rc{u>) = a + buj'^ , 



(6) 



we found that r is roughly equal to -1/2. From the defini- 
tion Eq.(0), we deduce that the critical Reynolds number 
Re is constant and is 



Re 



16. 



(7) 



C. Ripple wavelength (A) 

Using Eq.(||), we have studied the evolution of the 
mean ripple wavelength (A) from the center of the plate. 
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FIG. 6: Mean amplitude {A) as a function of the radius r. 
Different values of the water height h are illustrated. The 
solid vertical line gives the critical radius r^- One should note 
that Tc seems to be independent of h. 
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FIG. 7: Critical radius rc as a function of the rotation speed 
ci;. A fit using Eq.(^) is also illustrated. The power exponent 
is found to be equal to -1/2. 



to the sides of the container. Below the critical radius Tc, 
no modulation is observed and (A) is nearly zero. 

Within the unstable part of the interface, the measure- 
ment of the mean ripple wavelength is perturbated by the 
presence of "defects" . Indeed, when two ripples collapse 
a jump of (A) is recorded. We have noticed that on an 
image, the defects are almost situated at the same radius. 
Their role is to adjust the number of ripples to the shear 
conditions. In order to increase the accuracy of the mea- 
surements, we have decomposed the whole r-range into 
intervals which do not contain such "defects" . 

We have found that (A) grows linearly with the dis- 
tance form the center of the plate: 



{X) — (3 r for r > 



(8) 



This observation is consistent with the existence of spirals 
with constant arm numbers in the considered intervals. 
Moreover, we have noticed that the growth rate (3 does 
not depend on the height of water but decreases with the 
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rotation speed u, according to a power law with power 
exponent -1. Eventually, we see that the mean ripple 
wavelength obeys the scaling 

(A) - - for r > rc. (9) 

UJ 

In Figure ||, we have plotted /3 as a function of uj as well 
as a fit using Eq.(^). 
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FIG. 8: Log-log graph giving the growth rate /3 of the mean 
ripple wavelength (A) as a function of the rotation speed u. 
A fit using Eq.(S) is also illustrated. 



D. Azimuthal angle e 

Going deeper into the geometrical description of the 
patterns, we have measured the azimuthal angle e of the 
spiral arms. This parameter is defined at a given point by 
the angle between a ripple crest and the tangential line 
to the circle centered on the plate and passing by the 
considered point [see Fig||]. Practically, e is measured 
between two concentric circles with radii differing from 
Ar = 0.25 cm. 




FIG. 9: The azimuthal angle e is defined as the angle between 
the ripples crest (dashed line) and the tangential direction. 
The azimuthal angle e is measured between two concentric 
circles (solid lines) of radii r and r + Ar respectively (Ar = 
0.25 cm) 



An interpretation of the decrease of (3 can be deduced 
from the phase diagram illustrated in Fig.^. As the 
height of water is changed at fixed rotation speed, no 
major modification of the patterns is observed (nearly 
same wavelength). The consequence is that (3 does not 
depend on h. On the other way, we have seen that lj 
may cause the appearance of a second spiral. Such an 
occurrence results in a larger amount of spiral arms, i.e. 
in a larger number of ripples. Therefore, the mean ripple 
wavelength is lower when oj increases. The small values 
of (3 at large uj values clearly shows this fact. 

An important observation is that no well-defined tran- 
sition is observed in the curve giving (A) as a function of 
Lu. Indeed, the decrease of the mean ripple wavelength 
seems to be continuous. The way of appearance of the so- 
called second spiral is one reason explaining this result. 
At « 35 rpm, the second spiral is formed but has a 
smaller amplitude than the former spiral. The statistical 
weight of its wavelength is thus smaller than those of the 
former spiral. As the rotation speed increases, the sec- 
ond spiral amplitude growth. Eventually, the measured 
mean ripple wavelength corresponds to the number of 
spiral arms of both spirals, i.e. becomes at least twice 
those measured when only one spiral was observed. One 
should note that since (A) is divided by a factor larger 
than 2, we see that the hydrodynamics also tends to in- 
crease the number of arms of each spiral. This has been 
mentioned here-above as the origin of the "defects" in 
the spiral patterns. 



In order to determine the influence of the experimental 
parameters h, uj and r on the azimuthal angle e , we have 
proceeded as follows. At a given distance from the center 
of the container, i.e. a given r- value, and a given w-value, 
we measured e for different values of the height of water h. 
This process is repeated all over the container diameter. 
Then, the height of water h was fixed, and the procedure 
explained above was applied for different values of the 
rotation speed uj. 

We have noticed that the azimuthal angle e increases 
with the distance from the center of the container as well 
as with the height of water h and the rotation speed 
LD. Moreover, the measurements collapse when they are 
plotted as a function of the product hujr [see Figjl^l. It 
appears that the experimental data could be considered 
as following a power law as well as saturating exponen- 
tially. We have chosen the exponential law in order to 
avoid the azimuthal angle to grow to infinity. The fit of 
the data sets shows that 

/ hujr\ 

£ = £m - ai exp , (10) 

V "2 / 

where Em = 47° ± 8° is the asymptotical value of the 
azimuthal angle at high hur values, ai and a2 being 
fitting parameters. 

From Eq.([l0|), one can see that the azimuthal angle e 
meets the horizontal axis, says e =0, for the critical hur 
product value pc 

Pc = a2 In (^-^ ^ 3 IQ-^ m^ s^^. (11) 
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FIG. 10: Semi-log graph of the azimuthal angle e as a function 
of the rotation speed u> multiplied by the height of water h and 
the radial distance r. A fit using Eq.([lo|) is also illustrated. 



Since the azimuthal angle is zero at this critical value 
Pc, concentric circles must be observed for hujr = and 
spirals with azimuthal angles e < for values of hcur < 
Pc- However, we have never observed such patterns. 

An interpretation of this can be deduced from the con- 
sideration of the fluid momentum P — mv, where m is 
the mass of fluid, and v its linear speed. Considering that 
p = 1 10^ kg m~^ is the density of the fluid, and that the 
radius of the container R = 0.11 m, one gets 

P = pttR^ hujr, (12) 

what contains the hujr term. This thus means that the 
value of the azimuthal angle e is governed by the fluid 
momentum, according to Eq.(|l^). Combining Eq.(p^) 
and Eq.([l^) we are able to predict the fluid momentum 
corresponding to concentric circles. This critical momen- 
tum Pc reads 

Pc = /07ri?Vc = 0.011 kg m s-\ (13) 

For fluid momentum P < Pc, spirals with negative az- 
imuthal angles must be created. 

The fact that no concentric circles are observed comes 
from the lack of granular transport for such small fluid 
momentum. Indeed, ripples are only created when the 
grains are transported in the fluid, and this is only pos- 
sible if the fluid shear stress exceeds a certain threshold. 
Typically, ripples are observed only for fluid momentum 
P > 0.1 kg m s-^ 

IV. DISCUSSION 

The instability encountered in our experiment creates 
patterns very similar to those observed in the case of a 
fluid flowing over a disk at rest |l|, |, |, g]. The instabil- 
ity resulting in such patterns takes place in the Bodewadt 
boundary layer, i.e. in the boundary layer over the bot- 
tom disk. The patterns we observe are also similar to 



those observed in the Thomas' experiment In the 
latter experiment the ripples are created by an accelera- 
tion of the container. The acceleration rate Aw is some- 
what analogue to the rotation speed uj presented here. 
Hereafter, our experiment and the Thomas one will be 
referred as "granular" one, while 01, |,| will be de- 
fined as "fluid" experiments. 

Despite the similarity between the patterns, it is not 
obvious that the spiral patterns observed here are of the 
same nature as the fluid ones or as the Thomas ones. We 
are now attempting to make some comparison in order to 
differentiate the different kinds of instabilities. It is im- 
portant to notice that Bodewadt flows are very compli- 
cated to create experimentally. As a consequence, there 
is a poor amount of data. 

The major difference between the "granular experi- 
ments" and the fluid ones is the dependence of the mean 
wavelength of the patterns on the distance from the cen- 
ter [see Tabic ^. In Thomas and our experiments, the 
mean pattern wavelength increases with distance from 
the center while it is nearly constant in fluid experiments. 
However, the order of magnitude of (A) for all experi- 
ments is about 1 cm. 

The values of the azimuthal angle of the spirals ob- 
served in fluid experiments range between 13° and 17°, 
what is nearly constant compared to the angles measured 
by Thomas (between 20° to 70°) and the angles we mea- 
sured (from 4° at low fluid momentum to 45° at higher 
fluid momentum). The most important observation is 
that e increases with the radius in our experiment, while 
it decreases in 

Measurements of critical Reynolds numbers show that 
Re. decreases with the rotation speed in all experiments. 
However, the dependence is very slight in fluid experi- 
ments, compared to what is observed in granular ones. 
Moreover, Rc ~ 40 in fluid experiments and is only 
around 16 here. 



Scaling 


Parameters 


Here Thomas[9] Bodewadt [ 




(A) 

e 

Re. 


^ r/uj ^ r ^ cte 
^ hujr 7.-0-5 ^ 
« 16 -tj-"-i «40 



TABLE I: Behaviors of the physical parameters (A) , e and 
Re. for different experiments. 

From the comparison of the orders of magnitude and 
the behaviors, we can conclude that the instability ob- 
served here is not the signature of an instability in the 
boundary layer. What we observed here are ripples, sim- 
ilar to those observed in coastal areas, but in a rotating 
frame. Thomas pointed out the same conclusions. In 
both cases the spiral shape formed by the ripple crests 
may be due to the presence of centrifugal forces and other 
effects caused by the motion of the grains in the fluid. 
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However, due to the differences between the experimental 
procedures followed in and here, the created patterns 
do not behave in the same way. 

V. SUMMARY 

An experimental study of the instability rising on a 
sand-water interface has been proposed. The formation 
of the spiral patterns have been shown and have been 
described with the help of a phase diagram. 

The properties of the critical radius where the insta- 
bility forms have been studied. The critical radius Tc 
is found to decrease with the rotation speed of the con- 
tainer, according a square law, giving a critical Reynolds 
number equals to Re^ ~ 16. 

The mean ripples wavelength is proportional to the dis- 
tance from the center of the container and inversely pro- 
portional to the rotation speed of the container. Study- 
ing the azimuthal angles of the spirals, we have shown 



that this angle is governed by the momentum of the fluid 
P = mv, where m is the mass of fluid, and v its linear 
speed. Prom fits of experimental data, the azimuthal an- 
gle is found to saturates at around 45° for very high fluid 
momentum. A comparison between instabilities rising in 
the Bodewadt boundary layer, as well as a comparison 
with an experiment of ripple formation in an accelerated 
frame, allows us to conclude that the observed patterns 
are not the signature of a fluid instability. The instabil- 
ity of the sand-water interface we observe are similar to 
natural ripples created in a rotating frame. 
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